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It is shown that the mechanism of order out of disorder is at work in the antisymmetric pyrochlore 
antiferromagnet. Quantum as well as thermal fluctuations break the continuous degeneracy of the 
classical ground state manifold and reduce its symmetry to Z3 x Z2 . The role of anisotropic symmetric 
exchange is also investigated and we conclude that this discrete like ordering is robust with respect 
to these second order like interactions. The antisymmetric pyrochlore antiferromagnet is therefore 
expected to order at low temperatures, whatever the symmetry type of its interactions, in both the 
classical and semi classical limits. 
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I. INTRODUCTION 

The role of geometrical frustration in magnetic systems 
is at present one of the open question in the physics of 
strongly correlated systems since it can lead to novel low 
temperature behaviors. One of the most studied model is 
the nearest neighbor antiferromagnetic Heisenberg model 
on the pyrochlore lattice. For classical spins, the ground 
state is known to have macroscopic degeneracy, prevent- 
ing any magnetic ordering at T=0. In real systems, 
this degeneracy is often removed, partially or totally, 
by additional interactions such as dipolar interactions^, 
2nd neighbor exchange^, single ion anisotropy^,, magneto- 
elastic coupling^... This explains why pyrochlore com- 
pounds often order at low temperature, but their order- 
ing temperature is usually much smaller than the para- 
magnetic Curie temperature 8p since Tjq is not related 
to the exchange interaction J, but to the additional small 
interaction. Another type of process which can remove 
the degeneracy is order by disorder: as proposed first by 
Villain^ the classical degeneracy can be lifted by fiuctua- 
tions, if fluctuations around these classical ground states 
select one of these states. Both classical and quantum 
fluctuations are able to produce order by disorder, trough 
their contribution to entropy or zero point energy respec- 
tively. In this paper we study the effect of Dzyaloshinsky- 
Moriya interactions in the Heisenberg pyrochlore model. 
In an earlier paper—, we have shown that the T=0 K 
degeneracy is completly removed only for one particu- 
lar choice of the DM vectors. Here we show that if the 
DM vectors are such that the T=0 degeneracy is not 
removed, both classical and quantum fluctuations select 
6 equivalent ordered configurations, leading to a well de- 
fined magnetic structure. 



II. DMI AND MORIYA'S RULES FOR THE 
PYROCHLORE LATTICE 

Dzyaloshinsky has shown^ that, in crystals with no 
inversion center, the usual isotropic exchange J Si.Sj 
is not the only magnetic interaction and antisymmet- 




FIG. 1: The pyrochlore lattice. The interplay of the frustra- 
tion of the tetrahedral unit cell with the weak connectivity 
of the (corner sharing) tetrahedra provides peculiar magnetic 
properties to the pyrochlore antiferromagnet. 



ric exchange Dij.{Si x Sj) is allowed. In magnetic oxy- 
des, exchange interactions are usually attributed to su- 
perexchange mechanism, which involves virtual hopping 
between 2 neighboring magnetic ions. Taking into ac- 
count Coulomb repulsion and Pauli's principle, Anderson 
has explicited a microscopic mecanism which leads to 
isotropic super-exchange interactions. Later, Moriya has 
shownii-^ that inclusion of spin-orbit coupHng on the mag- 
netic ions in 1st and 2nd order leads to antisymmetric ex- 
change and anisotropic exchange respectively. Moriya's 
microscopic derivation of the DMI is only vaHd for insula- 
tors but other possible microscopic mecanism relevant for 
other materials were explicited, for instance in systems 
with RKKY interaction^ilS'. 

Whatever the microscopic origin of the DMI is, there 
are always symmetry constraints on the possible D vec- 
tors which may appear in the hamiltonian. Indeed, the 
hamiltonian must be invariant under the symmetry op- 
erations of the crystal and this will restrict the possible 





FIG. 2: D vectors for the DMI in the pyrochlore lattice. The 
convention is taken to fix the order for the cross products (al- 
ways D.Si X Sj with j > i). The two possible DMI are those 
obtained by varying the direction of the D vectors (D -^ — D). 
The case with the D as represented on the left is refered to 
as the "direct" case and the other case (on the right) is the 
"indirect" case. Once the DMI between two spins is fully spec- 
ified, the others DMI are also fixed and obtained by applying 
the different ^ rotations around the cube's diagonals which 
leave the tetrahedron invariant. The DMI in the rest of the 
lattice are also fixed and obtained by applying appropriate 
symmetry operations of the pyrochlore lattice. 





FIG. 3: One possible ground state in the case of indirect DMI. 
The ground state for the whole pyrochlore lattice is a q = 
structure so that only one tetrahedron is represented. Similar 
structures in the zx and yz planes are degenerate. 



coplanar phases, parametrized by one angle tp as follow. 



D vectors to those for which Dij.{Si x Sj) is invariant 
under these symmetry operations. This way of constrain- 
ing the D vectors has been given the name of Moriya's 
rules. Using these rules, we have shownii that two types 
of DM's interactions are compatible with the pyrochlore 
space group symmetry. These two patterns of D vectors 
have been given the name of "direct" and "indirect" case 
and are described in figure [2l 

As shown in our previous work, the classical ground 
state is always g = 0. It is non degenerate in the direct 
case and is described by the so called all-in all-out mag- 
netic phase; this case will not be discussed in this paper. 
From now on, we will focus on the second case of indirect 
DMI, where the ground state manifold is reduced com- 
pared to the case with no DMI, but is still continuously 
degenerate as shown in Ref. [ill. The role of anisotropic 
exchange will also be discussed in this paper. 



III. ROLE OF FLUCTUATIONS WITHIN THE 
GROUND STATE MANIFOLD 



From now on, we will focus on the second case of in- 
direct DMI, where the ground state manifold is g = 0, 
but is continuously degenerate. There are two kinds of 
ground states^ : (i) phases can be coplanar, belonging 
either to the xy, yz or zx plane and consists in perpen- 
dicular pairs of anticollinear spins (see figure [3]) which can 
rotate freely within these planes, (ii) There are also non 
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where the spins are labeled as in figure El and where (y? 
and 9 are not independant : 
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Consequently, it is possible to picture the classical 
ground state manifold by following the paths defined by 
the arrow of say, spin 1, in all allowed configurations (ac- 
tually, fixing any spin determines the configuration of the 
3 other spins). Such paths are depicted in figure HI 

The ground state is therefore a multiply arcwise con- 
nected one dimensional manifold, represented by the 
thick black fines on figure HI Six configurations play a 
particular role, since they are at the intersection of two 
lines, i.e they belong to both kinds of ground states de- 
fined above, where spins are coplanar, either in the xy, yz 
or zx plane (see Fig.[3l). In these phases, take for instance 
the xy plane, two spins are anti parallel (aligned along 
the (110) axis) and perpendicular to the two remaining 
spins (aligned along the (110) axis). These configura- 
tions, related by time reversal symmetry and rotations, 
are globally invariant under Z3 x Z2, and will be referred 
to this way from now on. 

It is of practical and theoretical interest to test wether 
fiuctuations operate on that manifold and reduce the de- 
generacy. Before getting into a quantitative analysis, it 
is worth pointing out qualitative arguments towards an 




A. Thermal fluctuations 

As emphasized in the previous paragraph we can guess 
that the pecuHar points of the configuration space around 
which thermal fluctuations could be eflicient are the 
above mentioned 6 points because they are the only one 
within the ground state manifold having two obvious "es- 
cape lines" e.g, two zero energy modes. It is therefore 
natural to deflne a parameter which is maximum and 
equal to one only when spin conflgurations are among 
the previous six. In each of these configurations, spins 
all belong to the same plane, either xy, yz or zx. When 
they belong to the xy plane, spins 1 and 4 are collinear as 
are spins 2 and 3, with the 1-4 pair being perpendicular 
to the 2-3 pair (see figure [31). The same occurs for the 
two other planes. Therefore, for each plane selection, a 
global index a can be defined to label a vector ef on each 
site of the lattice, so that the quantity 
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FIG. 4: Description of the one dimensional ground state man- 
ifold. Any move within the ground state can be described by 
the possible path followed by only one of the four spin as they 
are all parametrized by the same angle (see Eq.[T|). Here, only 
spin 1 is described. It can belong either to the xy, yz or zx 
plane, or be in a non coplanar phase. These non coplanar 
phases intersect coplanar phases at 6 points and the only way 
to go from one plane to another is to pass through these in- 
tersections, depicted by spheres. 



equals ±1 for the two time reversal related phases of each 
plane selection a = xy, yz, zx. If the spins belong to the 
xy plane, they point along the (1,-1,0) direction, thus 
^xy — ^(Ij ~1: 0) cind similarly for the two other planes. 
Defining then 
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order out of disorder mechanism. The main ingredient 
for fiuctuations to be efficient resides in the spectrum of 
the system : if two degenerate phases have different spec- 
tra and in particular, if one has softer modes, it is com- 
mon that this one is selected at low temperatures, where 
entropic effects and/or quantum fiuctuations may take 
place. The measure of the ground state manifold also 
plays an important role and should be sufficiently small. 
Schematically, if one phase is to be selected by fiuctua- 
tions but is "drowned" in a huge ground state manifold, 
the system cannot statistically pick up that phase and 
fiuctuations will therefore not operate. In the present 
case, as it is shown hereafter, both preceding conditions 
are fulfilled : there are obviously peculiar points expected 
to have spectra with softer modes (zero modes actually) 
and the ground state manifold, because it is one dimen- 
sional, has a null measure in configuration space and al- 
lows order out of disorder to be efficient, both at thermal 
and quantum levels. 



allows for characterizing any spin configuration and its 
proximity to one of the six peculiar phases because M = 
1 for these configurations and these configurations only. 

In order to test wether thermal fiuctuations may en- 
tropically select peculiar phases of the ground state mani- 
fold, we have performed classical monte carlo simulations 
of finite size lattices with periodic boundary conditions. 
The clusters we investigated had 32, 500 and 2048 sites. 
For each simulations, we used a single fiip metropolis al- 
gorithm with a modified update so that the acceptance 
rate stayed roughly above 40%, associated to a local ro- 
tation of each spin around its local molecular field. Dur- 
ing the simulation, the autocorrelation time is calculated 
on the fiy in order to adapt the number of monte carlo 
steps between two measures and ensure that measures 
are uncorrelated. Results are reported in Fig [5] and [6] 
for respectively D / J = 0.1 and D / J = 1. The measure 
of (to) as a function of the temperature clearly indicates 
that this parameter is an order parameter, although we 
cannot discuss the order of the transition so far. 

Both cases support that at low temperature, the sys- 
tem orders in one of the 6 previously described states as 
was expected from qualitative arguments. Moreover, it 
can be seen that Tc is of the order of D, which confirms 
that DMI are responsible for this ordering. 
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FIG. 5: Thermal average of m (see Eq. [3]) as a function of 
temperature for three lattice sizes, N = 32, 500 and 2048, 
and for D/J = 0.1. 
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FIG. 6: Thermal average of M (see Eq. [3|) as a function of 
temperature for three lattice sizes. A*' — 32, 500 and 2048, 
and for D/J = 1.0. 



B. Quantum fluctuations : semiclassical approach 

The same investigation has been done for quantum 
fluctuations. Contrary to the classical case, it is not nec- 
essary to guess which points will play a particular role. 
This is because it is possible to compute at zero temper- 
ature the energy renormalized by quantum fluctuations 
for all phases of the ground state manifold. To do so, 
we have performed a linear spin wave expansion of the 
Hamiltonian around each state of the ground state mani- 
fold, e.g, around each state parametrized by the angle ip. 
The resulting bosonic hamiltonian is then diagonalized 
which allows to obtain the normal modes and calculate 
the renormalized energy and magnetization by quantum 
fluctuations. 

The first result is that the Z3 x Z2 phases are unam- 
biguously selected by quantum fiuctuations, as shown in 
Fig. [71 For planar phases, they correspond to a = or 



FIG. 7: Renormalized energy per site as a function of the 
angle a, for planar phases [all planes are equivalent], with 
D/J = 1. a = and a = tt correspond to one time reversal 
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FIG. 8: Renormalized magnetization for S = 1/2 as a func- 
tion of the angle a, for planar phases, with D/J = 1. a = 
and a = n correspond to one time reversal related pair of 
states belonging to the Z3 x Z2 phases. At these two points, 
fluctuations have softer spectra, hence local magnetization is 
much more renormalized. 



a = TV, a being the angle of the spins with one of the 
6 peculiar phases (2 for each plane). The same calcu- 
lation was performed for conic phases (indexed by (p), 
confirming that the 6 phases are selected by quantum 
fluctuations. 

This selection is also illustrated by the computation 
of the renormalized magnetization, as shown in Fig. [8l 
As conflrmed by the calculation of the energy, there are 
softer modes near the 6 selected phases. Consequently, 
transverse fluctuations are much stronger giving a larger 
renormalization of magnetization around these phases. 

Wether the magnetization, which is zero for D/J = 0, 
could be stabilized by taking DM's interactions into ac- 
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FIG. 9: Renormalized magnetization of one of the symmetry 
related Z3 XZ2 phases, as a function oi D/J for S — 1/2, 1, 3/2 
and 2. For S* > 1, a DM interaction of 1 percent is sufficient 
to stabilize a Neel like long range order. From that point of 
view, DM interactions are very efficient in destroying the clas- 
sical spin liquid behavior and driving this model to a classical 
behavior. 



dipolar interactions^, 



K- 
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where indices i,j are related to the fee Bravais lattice 
while indices m, n correspond to one of the four sites of 
the tetrahedral unit cell. The set of unit vectors n" is 
(1,0,0), (0,1,0), (0,0,1). We allow the sign of AT to vary as 
there is no obvious reason why one peculiar type, K > 
or K < 0, should prevail. The full hamiltonian we are 
left with reads now, 
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where J^ (Rij^) are the coefficients of the coupling ma- 
trix J {q) which contains the isotropic symmetric ex- 
change —JSi.Sj, the antisymmetric exchange Dij.Si x 
Sj and the anisotropic symmetric exchange K (see 
Eq. SJ. Stabilized low temperature phases are investi- 
gated through a mean field analysis. The hamiltonian is 
rewritten in reciprocal space using the following trans- 
formations, 



count is also of interest because it quantifies how eas- 
ily this model becomes "classical". We have therefore 
computed for one of the six selected phases, the value of 
the quantum renormalized magnetization as a function of 
D/J. Results are reported in Fig. [9l It is worth noting 
that apart from the case S — 1/2, a very small value of 
D/J {D/J > 0.01) stabilizes the magnetic ground state 
of the pyrochlore antiferromagnet. Furthermore, we can 
quantify how efficient DM's interactions are by fitting the 
gain in energy as a function of D/J. For small values of 
D/J, quantum fluctuations contribute to the renormal- 
ization of the energy as \JD / J . This is an important 
point as discussed in the next section. 



IV. ROLE OF ANISOTROPIC SYMMETRIC 
LIKE EXCHANGE 



So far, only DM's interactions have been taken into 
account. Whatever the type of fluctuations, thermal or 
quantum, the order out of disorder mechanism is effi- 
cient and reduces the continuous degeneracy to a discrete 
one. But is this selection robust vis a vis next order (in 
spin orbit coupling) interactions present in the Hamil- 
tonian, the anisotropic symmetric exchange? Unfortu- 
nately, we must lose generality to address this question 
because many kind of symmetry allowed interactions are 
possible. It is nevertheless possible to argue quantita- 
tively on that matter. 

In order to be specific, we choose to consider inter- 
actions that keep all symmetries of the pure pyrochlore 
lattice. They correspond to nearest neighbour truncated 
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where N is the number of Bravais lattice points. The 
resulting interaction matrix J{q) is a 12 x 12 non diag- 
onal hermitian matrix. Hence, to completely diagonalize 
J{q) one must transform the q-dependent variables, S™, 
to normal mode variables. In component form, the nor- 
mal mode transformation is given by 



p— 1 7— 1 
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where the indices (p, 7) label the normal modes (12 for 
Heisenberg spins), and {(/>q'''} are the amplitudes of these 
normal modes. U(q) is the unitary matrix that diagonal- 
izes J{q) with eigenvalues A(q). Finally, the mean-field 
free energy to quadratic order in the normal modes reads, 
up to an irrelevant constanlii^. 
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where J-{T) is the mean-field free energy per unit cell, T 
is the temperature in units of fce , and n = 3 for Heisen- 
berg spins. Therefore, the mean-field low temperature 
phase is defined by the corresponding wave vector Qord as- 
sociated with the extremal eigenvalue maxp,^_q (A^(q)). 
When iiT = 0, it was shown in preceding sections that 
Qord = and that the ground state manifold is con- 
tinuously degenerate. Including non zero K, whatever 



its sign, the degeneracy is lifted already at the mean 
field level. The new ground states are, for both posi- 
tive and negative K, q — slightly distorted version of 
the Z3 X Z2 states. The energy decrease, for small K, is 
always quadratic, AE ex K^. 

At this point, we can discuss wether anisotropic sym- 
metric interactions may interfere with quantum fiuctua- 
tions induced by DM's interactions. The first remark is 
that at the mean field level, K breaks the degeneracy but 
select phases which are continuous deformation of quan- 
tum selected phase. Therefore, we expect that taking 
into account anisotropic exchange would not change the 
behavior of this system. The second point is related to 
how strong K is efficient in selecting a particular phase. 
Suppose that the mean field selected phase is not one 
of the Z3 X Z2 phases or a slightly distorted version of 
those. Because the energy gain is quadratic in K, it ap- 
pears that quantum fiuctuations are much more efficient, 
as they induce an energy gain proportional to Vd (see 
section llllBp . This clearly shows that even if anisotropic 
interactions are taken into account, we expect a lifting 
of degeneracy corresponding to the Z3 x Z2 phases or 
slightly distorted versions of those. The only mechanism 
which could work against such scenario would be that 
quantum fiuctuations around K selected phases be larger 
than ^/D. Because K (x D'^, this would mean that one 
should expect the energy gain induced by K-like interac- 
tions to behave like vK, which is highly unprobable. 



romagnet. When the symmetry type of the allowed DMI 
leads to a one dimensional degenerate ground state, it is 
shown that these two types of fiuctuations are efficient 
in reducing the degeneracy and drives the model to the 
same ordered ground state with a discrete global degen- 
eracy of the Z3 X Z2 type. Including higher order terms in 
the hamiltonian like the anisotropic symmetric exchange 
is likely to leave thermal and quantal fiuctuations dom- 
inate the low temperature behavior of this system. It is 
therefore concluded that whatever the type of symme- 
try allowed DMI, direct or indirect, they should drive 
the pyrochlore Heisenberg antiferromagnet to a low tem- 
perature ordered phase. When the model is essentially 
classical, the critical temperature is of order D. If quan- 
tum fiuctuations are present, the value of the ordering 
temperature results from a balance between the strength 
of the DMI and the magnitude of the fiuctuations. As 
we have shown that for moderately large DMI the local 
magnetization is quite stable, the value of the ordering 
temperature should not be very different from the one of 
the classical case. Finally, it is worth pointing out that 
most of studied pyrochlore compounds are Rare Earth 
oxides. It is therefore the next natural step to investigate 
the mechanism of DMI within 5f alloys and determine 
whether those interactions could be in some compounds 
large enough to drive the low temperature behavior of 
these frustrated systems. 



V. CONCLUSION 

This works investigates the role of thermal and quan- 
tal fiuctuations on the antisymmetric pyrochlore antifer- 
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